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Abstract. The barrier heights (including zero-point effects) for H + CHoO — * CH 3 O 
and CH 3 O —* CH 2 OH have been computed using complete active space self consis- 
tent field (CASSCF)/ gradient calculations to define the stationary point geome- 
tries and harmonic frequencies and internally contracted configuration-interaction 
(CCI) to refine the energetics. The computed barrier heights are 5.6 kcal/mol and 
30.1 kcal/mol, respectively. The former barrier height compares favorably to an 
experimental activation energy of 5.2 kcal/mol. 
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The CH 3 0 radical is important in combustion and atmospheric chemistry [ 1 ]. In 
addition, recent experiments have studied individual rovibrational states of CH 3 O 
above the H-CH 2 0 dissociation threshold by stimulated emission pumping [3]. The 
interpretation of these experiments depends critically on the barriers to decom- 
position of CH 3 0 to H + CH 2 0 and rearrangment to CH 2 OH. Previous studies 
of these barrier heights include the work of Saebp, Radom, and Schaefer [2], who 
used M 0 ller-Plesset perturbation theory through third order with small basis sets 
up through 6-31G**. More recently, Page, Lin, He, and Choudhury [4] reported 
somewhat more accurate calculations using multireference configuration-interaction 
(MRCI) with polarized triple zeta basis sets. In the present paper more accurate 
calculations are reported which make use of larger basis sets and extensive Cl. 

Two different basis sets were used in this work. For the complete active space 
self consistent field (CASSCF) gradient calculations the polarized double zeta set 
of Dunning and Hay [5] was used. The basis set for C and 0 is a (9s5p)/[3s2p] 
basis augmented by a single set of 3d functions with exponents of 0.75 and 0.85 
for C and 0, respectively. The H basis is (4s)/[2s] augmented with a single set of 
2 p functions with exponent 1 . 00 . The basis set used in the Cl calculations is the 
Dunning correlation consistent triple zeta double polarization basis set [ 6 ]. This 
basis is [ 4 s 3 p 2 dlf] for C and 0 and [ 3 s 2 pld] for H and is described in detail in Ref. 
6 . 

The calculations were carried out in C, symmetry for a wavefunction of "A 
symmetry. From the Hessian matrix in the diagonal representation, the station- 
ary points obtained in C, symmetry are the minimum for CH 2 0 -I- H, the saddle 
point for CH 2 0 + H -» CH 3 0, the CH 3 0 minimum, and the saddle point for 
CH 3 0 —► CH 2 OH. For CH 3 0 the 2 A' state corresponds to one component of the 
2 E state in C 3 U symmetry. According to the calculations of Saebp et al. [2], the 
Jahn-Teller stabilization energy in this system is 0.56 kcal/mol and the Jahn-Teller 
splitting is 0.12 kcal/mol, with the 2 A' component lower. CH 2 OH is found to have 
no symmetry and the use of a mirror plane of symmetry would not be appropriate. 
However, the separation between CH 3 O and CH 2 OH has been accurately computed 
by Bauschlicher, Langhoff, and Walch [7] and further calculations for CH 2 OH were 
not carried out here. 
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The CASSCF calculations consisted of five electrons and five orbitals. The active 
electrons for CH 3 O included the CO bond pair, one CH bond pair (the CH bond in 
the molecular plane), and the a' 0 2 p like orbital. These orbitals correspond to the 
CO <7 and 7 r bonds plus a H Is orbital for CH 2 0 + H. The remaining electrons, which 
are inactive, include the other two CH bond pairs, the a" O 2p like orbital and the 
0 Is, 0 2s, and C Is like orbitals. In generating the set of reference configurations 
for the subsequent internally contracted Cl (CCI) calculations, no more than two 
electrons were permitted in the weakly occupied CASSCF orbitals. All but the 0 Is 
and C Is electrons were correlated in the CCI calculation. A multi-reference analog 
of Davidson’s correction [ 8 ] was added to the CCI energies. 

The CASSCF/gradient calculations used the SIRIUS/ABACUS system of pro- 
grams [9], while the CCI calculations were carried out with MOLPRO [ 10 , 11 ]. The 
calculations were carried out on the NASA Ames Cray Y-MP. 

Table I gives the harmonic vibrational frequencies, rotational constants, and ge- 
ometric parameters for the stationary points considered in this work. Here H a is 
the H of the spectator CH bonds, while H*» is the hydrogen which transfers during 
the reaction. The geometries and frequencies in Table I are quite similar to those 
obtained by Page et al. [ 4 ]. This is expected since a CASSCF wavefunction was 
used in Ref. 4 . However, in Ref. 4 the CO cr bond was not correlated. For CH 3 0 
this has the expected effect that the CO bond is elongated by ^ 0.03 A compared to 
the calculations in Ref. 4 . At the H-CH 2 0 saddle point the CO bond is elongated 
by ^ 0.02 A and the C-H& bond is shorter by ^ 0.05 A as compared to Ref 4. 
There are larger differences between this work and the results of Saebp et al. Here 
the geometry for CH 3 O is similar except for the expected elongation of the C-H^ 
and CO bonds. However, there are significant differences in the geometry of the 
H-CH 2 0 saddle point. In particular, the C-Hj, bond length is shorter by « 0.18 A 
. This result is consistent with the expectation that saddle point geometries will be 
more sensitive to electron correlation effects than will stationary points. 

Table II gives the computed energies and zero-point corrections. The energies 
are from the CCI calculations, while the zero-point corrections are derived from the 
CASSCF harmonic frequencies. From Table II it is seen that the relative energies are 
only very slightly changed (« 0.1 kcal/mol) in going from the [ 4 s 3 p 2 dlf/ 3 s 2 pld] 


3 



basis set to the [5s4p3d2f/4s3p2d] basis set. Table III compares the computed 
barrier heights obtained in the present calculations to those obtained by Saebp et 
al. [2] and by Page et ad. [4]. Here it is seen that the barriers obtained by Saebp 
et al. are too large by « 6 kcal/mol., while the barrier for H + CH2O — + CH3O 
obtained by Page et ad. is 2.3 kcal/mol larger than in the present calculations. Our 
computed barrier height is 5.6 kcal/mol for H + CH2O — * CH3O and 5.2 kcal/mol 
for D + CH 2 0 — ► CDH 2 0. This result is in good agreement with experimental 
estimated activation energies of 5.2 kcal/mol for H + CH2O and 3.9 kcal/mol for 
D + H 2 CO [12]. 

It is concluded that the present calculations give energetics for the CH 3 0 sys- 
tem which are of chemical accuracy ( in this case the estimated error bars are 
± 1 kcal/mol). This is in contrast to results obtained by Saebp et al. using Mpller- 
Plesset perturbation theory, where the errors in the computed barrier heights are 
% 6 kcal/mol. 
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Table I. Computed stationary point harmonic frequencies (cm *), rotational 
stants (cm -1 ), bond lengths (A), and bond angles (degrees). 
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Table II. Computed energies and zero-point corrections. 
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Table III. Comparison of computed barrier heights (kcal/mol). 



Saebp et al. a 

Page et al. 6 

PW C 

Exp. 

H + CH 2 0 CH 3 0 

12.5 

8.0 

5.6 

5.2 d ,6.2 

D + CH 2 0 -*• CDH 2 0 



5.2 

3.9^ ,4.8 

ch 3 o -> CH 2 OH 

36.0 


30.1 



“ Reference 2. 
b Reference 4. 
c Present work. 

d Estimated activation energy from a room temperature rate constant for D + 
H 2 CO, an assumed preexponential factor of 2 x 10 13 cm 3 mol 1 sec 1 , and a zero- 
point energy difference between D + H 2 CO and H + H 2 CO of ~ 1.3 kcal/mol. Ref. 
12 . 

e Estimated threshold energy from an RRKM analysis. Ref. 12. 
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